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Abstract 

Optical potentials for elastic p-d scattering and the coupled processes p+ 3 He 
— > p+ 3 He and p+ 3 He — > d+d are derived in the Faddeev-Yakubovsky frame- 
work with special emphasis on leading order terms, which are expected to be 
valid at intermediate energies. In addition, equations for the fragmentations 
3 He(p,ppp)n and 3 He(p,pp)d are derived within the same framework. Again 
leading order terms for intermediate energies are considered. 
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I. INTRODUCTION 



Optical potentials have a long tradition in describing scattering of protons and 
neutrons from composite targets. At intermediate energies their theoretical formulation is 
often based on approaches like the Watson 0, the KMT formulation or the spectator 
expansion 0-§J. The fundamental idea from which those approaches start is the grouping 
of the scattering process for a nucleon hitting a nuclear target into rescatterings of various 
orders. In lowest order the projectile interacts with one target nucleon, which has a mo- 
mentum distribution according to a mean field generated by the residual nucleus. In second 
order the struck target nucleon rescatters from a second target nucleon, which therefore par- 
ticipates actively in the scattering process and is thus no longer 'hidden' in the mean field 
of the target. This more and more complex scenario can be formulated in a general fashion 
in the so called spectator expansion, but in practice only the lowest order processes have 
been numerically realized. In the most elaborate case this leads to the so called full-folding 
model |0-§[. Here, for constructing the optical potential the fully-off- shell nucleon-nucleon 
(NN) t-matrix is convoluted with the single nucleon density matrix of the target. Anti- 
symmetrization is kept among the target nucleons, but the projectile nucleon is treated as 
distinguishable and only the antisymmetrization between the active nucleons is kept. This 
procedure of treating the antisymmetrization is clearly approximative and there might be 
some flaws in the actual realization as pointed out in the Appendix A. However, this picture 
turned out to be quite successful in describing proton and neutron scattering from light ( 16 0) 
to heavy ( 208 Pb) nuclei at intermediate energies using a NN t-matrix derived from a realistic 
NN force and single nucleon density matrices resulting from nuclear structure calculations 
0-0. Since the structure part faces solving the full A-body problem, approximations have 
to be made. In addition, it is not the same NN force used in the structure part as the one 
acting between projectile and target nucleon. This can be considered as an inconsistency. 
The same theoretical approach has been used to describe elastic scattering of protons from 
a light nucleus like 3 He at intermediate energies P,|10 . 



At present only for very light nuclei there is a chance to eliminate this inconsistency 
concerning the nuclear structure part and treat the scattering process as well as the structure 
on an equal footing, namely based on the same realistic NN forces. It is the aim of the present 
study to derive rigorously the optical potential and transition potentials for scattering of 
nucleons from the two lightest nuclei, the deuteron and 3 He within the framework of the 
Faddeev-Yakubovsky equations. In this framework, antisymmetry between all nucleons is 
treated correctly. To exhibit differences as well as similarities to the spectator expansion, 
we give a brief sketch of the latter scheme in Appendix A. 

In Section II we consider the proton-deuteron (p-d) system and in Section III the 
more complicated p- 3 He system. In the latter case it appears natural to generalize the optical 
potential for elastic p- 3 He scattering to a potential matrix corresponding to a coupled set 
of equations describing scattering in both two-body fragmentation channels, p- 3 He and d-d. 
An additional interest lies in the study of the various breakup processes in p- 3 He scattering 
into 3 fragments, proton-proton-deuteron (ppd), and into 4 fragments, pppn. For these 
reactions a wealth of experimental information is available [0,0] • However, these data are 
conventionally analyzed in PWIA or DWIA only, which does not correspond to a microscopic 
four particle theory based on NN forces. Three- and four-body fragmentation processes 
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will be described in Section IV. Throughout the whole article we only consider the strong 
interaction when describing protons and neglect Coulomb forces. In a practical calculation 
the Coulomb interaction between the projectile proton and the charge distribution of the 
nucleus can be included in a straightforward manner at least in the elastic channel []10yi3l . 

In all cases to be discussed, we emphasize intermediate energies where the optical 
potential expressions and the breakup amplitudes can be expected to simplify. Presently, 
in the four-body system these simpler expressions appear to be the only ones which can be 
realized in practical applications. Nevertheless, they carry important information on the 
reaction process and can serve as testing ground for applying the free NN force in such 
reactions and as a means to study properties of realistic 3 He wave functions. Appendix B 
contains the leading terms of the transition potentials between the p- 3 He and the d-d chan- 
nels. 

The four-body scattering problem has been formulated in some detail in the con- 
text of the Alt-Grassberger-Sandhas (AGS) equations [rj,TJ|. This formulation is basically 



equivalent to the Faddeev-Yakubovsky one which we are using. In practice however, and 
especially in the context of deriving optical potentials for p- 3 He reactions, which is our aim 



here, it is not necessarily convenient to start from the equations presented in Ref. |L4|. In 
addition, we want to keep our formulation totally independent of assumptions like finite 
rank approximations to the underlying force. For these reasons we find it adequate to set 
up the coupled four-body equations our way by regarding the asymptotic behavior of the 
Faddeev-Yakubovsky equations in configuration space and to introduce at the earliest stage 
the identity of the nucleons. For the best of our knowledge, this sort of derivation has not 
been displayed before for scattering processes. For four-nucleon bound states however, this 



kind of formulation has already been used [16]. We conclude in Section V. 



II. ELASTIC PROTON-DEUTERON SCATTERING 

It is well known that the operator U for elastic p-d scattering obeys the AGS 



equations [|l4l ,|T5|[ 



U = PG 1 + PtG U. (2.1) 

Here Go is the free three-nucleon propagator, t the off-shell NN t-operator and P a sum of a 
cyclic and an anticyclic permutation of 3 particles. This operator equation has to be applied 
onto the initial channel state 



|$ qo >= \4> d > |q > (2.2) 

where \<pd) is the deuteron state and |qo) the momentum eigenstate of relative motion of 
the projectile nucleon with respect to the deuteron. Throughout this Section we use the 
standard Jacobi momenta p for a two-body subsystem and q for the 'spectator' particle. In 
addition we work in the isospin formalism and treat the nucleons as identical. 

The motion of the three nucleons in intermediate states is now separated into two 
parts, namely one, where a pair forms a deuteron, and another one where this pair is in 
the two-body continuum. Denoting by H the operator for the kinetic energy of the three 
nucleons and by V the potential operator for the NN force, we obtain 
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with 



tG = V{E + ie-H Q -V)- 1 = VG b + VG c (2.3) 



G b = ]</>*) „,. \ s-3 < H (2-4) 

E + is - E d - 



and 



G c = [ d*p\<p p >(+) — * < (+)| . (2.5) 

J E + ie - 2 -r-q 2 

m Am ^ 

The deuteron binding energy is given by Ed. The set of two-body eigenstates \<p d > and 
\(j)p >^ span the two-body Hilbert space and have been used in Eq. ( |2.4| ) and Eq. ( |2.5j ) to 



decompose the resolvent operator of Eq. ( |2.3|) . In Eqs. (2A) and (|2.5|) the denominator is 



still an operator acting on the spectator motion through the kinetic energy operator -£^q 2 - 
The driving term in Eq. (|2.1| ) when applied onto the state |$ qo ) given in Eq. (|2.2|) 



can also be written as 

PG l \% o >=PV\% o >, (2.6) 
where V is acting only within the deuteron. 

When inserting the decomposition given in Eq. ( |2.3| ) into Eq. Q2.1| ), we obtain 

U\% >= PV\^ > +PVG b U\% > +PVG C U\% > ■ (2.7) 

We now define 

U c \% >= PV\^ > +PVG C U C \% >, (2.8) 
and consequently can rewrite Eq. Q2.71) as 

U\% >= f/ c |$ qo >+J d 3 qU c \% > O (q) < ^\U\% > . (2.9) 



Here Eq. Q2.4|) has been written out explicitly as 

G b = fd 3 q\4> d > |q > — — g—j < <f> d \ < q| = /d 3 g|$ q > g (q) < $ q | . (2.10) 

J e + is — t, d — ^<r J 



With this preparation Eq. (|2.9| ) can be cast into the closed form of an integral equation by 
projecting onto the state ($ q /| 

< $ q '|C/|$ qo > = < ^l^cl^qo > +/ A < $ q /|C/ c |$ q > (/ (q) < $ql^l$qo > (2-H) 

This is the desired integral equation for the elastic p-d scattering amplitude. After a partial 
wave decomposition it becomes one-dimensional. The driving term 

<%>\U c \%>=V <ll(l (2.12) 
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can be interpreted as optical potential Vq/ q . The question of interest is, whether at inter- 
mediate energies (say above 100 MeV) one can expect that the integral equation, Eq. Q2.8|) , 
can be solved by iteration and moreover whether the very first few terms will be sufficient. 
Assuming this to be correct we expand 

< $q'|£/ C |$q >~< $q'|PV|$ q > + < %' \ PVG C PV | $ q > +■ ■ ■ (2.13) 

In order to proceed further, we need to be more specific and identify a two-body subsystem, 
e.g. we choose the pair (23) to be the deuteron. Then the potential V becomes V = V23. 
The permutation operator P is given by 

P = P12P23 + P13P23, (2.14) 

and the channel state by 

|$q >= |$q >1= \<f>d >23 |q >1- (2.15) 

The first term of Eq. ( |2.13| ) is now explicitly given as 

< $ q /|PV|$ q > = l< $ q '|1/ 3 l|$q >2 + 1 < $q'|K2|$q >3 , (2.16) 

where 

|$ q >2 = PwPnl&n >i= \<f>d >3i |q >2 (2.17) 

|$q >3 = As^sl^q >1= \4>d >12 |q >3 • 

With this, one easily finds 

< %<\PV\% > = 1 < $ q ^ 13 (l - Pl 3 )|$q >2 • (2.18) 

The potential V13, the NN force, should only contribute in two-body states which obey the 
Pauli principle. Therefore, it has to be represented by antisymmetrized states and thus 

V 13 (l - P 13 ) = 2V 13 . (2.19) 

It follows that 

< $q'|PV|$ q >= 2 x < $ q ,|T4 3 |$q >2 • (2.20) 

This is the well known single particle exchange process where the deuterons in the bracket 
states are formed out of nucleons (23) and (13), respectively. If the spectator momenta q 
and q' are larger than the typical momenta inside the deuteron this matrix element is small 
and will not be the leading contribution to the optical potential. 

Let us now consider the second term in Eq. ( |2.13| ). It is convenient to introduce a 
t-operator t c defined as 

VG C = t c G . (2.21) 
From its very definition it is related to the full t-operator by 
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t c = t - V\(j) d > g < <p d \G 1 (2-22) 

and obeys the Lippmann Schwinger equation 

t c = V-V\<j>i><4> d \+t c G Q V. (2.23) 

With the help of Eq. (|2.22| ) and Eq. ( p.l4j ) the second term in Eq. (|2.13| ) can be expressed 
as 

<^\PVG C PV\%> = < $ q ,|Pi c G PV|$ q > (2.24) 

= < $ q ^ c (31)(l - P 13 )G V 23 \% > 
+ <$ q ,|t c (21)(l-P 12 )G V 2 3| < f q >- 

If the momentum q would be on shell, i.e. |q| = |qo|, then we would have GoV^l^q) = 
|$ q ) and we would end up with a simple form, namely that nucleon 1 interacts with the 
constituents of the deuteron (nucleons 2 and 3) through t c . However, the expression is more 
complicated. The free propagator is given as 

Go = = \, - G (E -E d - ^q 2 ) 1 . (2.25) 

Therefore, 

G V 23 \^ >= \% > -CWl$q >, (2.26) 

so that the second term in Eq. (|2.13|) is given by 

<$ q ,|PKG c PV|$ q > = 
<$ q ,|t c (13)(l-P 13 )(l-G ^ 1 )l« , q> +<$ q '|t c (12)(l-P 12 )(l-Go^ 1 )l« , q> • (2-27) 

This expression can be reduced by replacing (1 — P^) by 2 and noticing that because of the 
antisymmetry of the deuteron states the two matrix elements are the same. This resulting 
form for the optical potential has some similarity to the often used 'tp' form as discussed in 
Appendix A. It will be interesting to perform numerical studies in order to see quantitatively 
the validity of the truncation in Eq. ( p. 13 ) for the optical potential and to compare this 
expression with the more standard 'tp' form. Such investigations are planned. 



III. TWO-BODY FRAGMENTATIONS IN PROTON- 3 HE SCATTERING 

Four-nucleon scattering has not yet been numerically mastered in a rigorous and 
general manner similar to the scattering of three nucleons upon each other ||17|| . However, 
at intermediate energies there is a chance, that specific approximations can be theoretically 
justified and systematically controlled in a numerical realization. This can open interesting 
insight into the reaction itself, the use of the free (not modified through the presence of 
the nuclear medium) NN interaction in such processes as well as into the properties of the 
three-nucleon bound state. Our approach is a rigorous framework and has been successfully 
used for the a-particle bound state |18|J19[| . 
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As is well known the fully antisymmetric scattering state for p- 3 He scattering, 
obeys the homogeneous integral equation 



¥« = G 5><,¥« (3.1) 



where Go is the free four-nucleon (4N) propagator and Vij represents the NN force between 
particles i and j. As a first step one usually decomposes into 6 Faddeev components 

*«=5>y (3.2) 

with 

% = GoVijVW = GoVij ]T ipki- (3.3) 

kl 

When introducing the off-shell NN t-matrix t{j, as already defined in Section II, one arrives 
at 6 coupled Faddeev equations 

ipij = GqUj ^ki- (3.4) 

kl^ij 

In the following we immediately make use of the identity of the four particles. 
This simplifies the notation significantly and it appears not to have been outlaid before in 
the literature in this particular form. It also serves to clearly define our notation. Without 
loss of generality we consider the scattering of particle 4 from a target consisting of the 
subcluster built from particles 1, 2 and 3. 
Let us consider the Faddeev component 

V>12 = G t 12 (^23 + fal + ^24 + V>14 + ^u), (3.5) 

which is then split into 3 Yakubovsky components, namely 

V>i = £0^12(^23 + fai), (3.6) 



Pu^i = -^0^12(^24 + i/>u), (3.7) 

and 

^2 = Gofaipu- (3.8) 

To arrive at Eq. (|3.7| ) we used the antisymmetry of \&( + ) and the definition of the Faddeev 
components as given in Eq. ( |3.3| ). This gives 

^12 = (1-P 3 4)^i + ^2- (3.9) 
With the same reasoning one has 

^23+^31=^12, (3.10) 
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where the permutation operator P is given in Eq. (|2.14j) . Thus Eq. (|3.6| ) can be written as 
V>i = G t 12 P^j 12 = G t 12 P((l - P 34 )^i + fa)- (3-11) 
The decisive step to describe p- 3 He scattering is to sum up all pair forces within the 3-body 



subcluster of particles 1,2 and 3. In order to achieve this, we rewrite Eq. Q3.ll 



as 



(1 - GotnP)^ = G t 12 P(-P u ^ + V 2 ). (3.12) 

Now the left hand side alone has a nontrivial solution, which is related to the incoming 
channel state via 

(1 - G t 12 P)$ F = (3.13) 

with 

$ F = |0 F (123))| U >4. (3.14) 

Here \(p F ) denotes the Faddeev component to the target state, namely the 3N bound state 
10), 

|0>=(1+P)|^>. (3.15) 

The momentum eigenstate of the projectile, |u > 4 , is described by an appropriate Jacobi 
momentum u. With these definitions we can write Eq. ((3.12 ) as 



fa = $ F + (1 - GotuPy'Got^Pi-Pu^ + ^2). (3.16) 

Defining now 

(1 - G Q t 12 P)- l G G t l2 P = G TP, (3.17) 
which is equivalent to introducing the 3-body operator 

T = t 12 + t 12 PG T, (3.18) 

one obtains 

= ® F + G TP(-P 3 4^i + V> 2 ). (3.19) 

The three-body operator T given in Eq. ( |3.18| ) is defined by an off-shell Faddeev equation 
for the three-body subsystem composed of particles 1,2, and 3, and depends parametrically 
on the kinetic energy of the fourth particle. The above Eq. fl3.19p is the first one of the 
Yakubovsky-equations. 

Let us now consider the second Yakubovsky component from Eq. (|3.8|) . With the 
help of the permutation operator 

P = PisPm (3.20) 

we can write 
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^2 = GotuPipu = G t 12 P((l - Pm)4>i + V>2)- 



(3.21) 



Next, we sum up the pair forces in the two noninteracting two-body subsystems (12) and 
(34), representing two deuterons, to infinite order. To do this, we rewrite Eq. ( |3.21| ) as 

(1 - GotP)tp 2 = G tP(l - P 34 )Vi (3.22) 

and solve for ip 2 . In the two-deuteron channels there are no ingoing waves, therefore the 
nontrivial solution to the left hand side alone should not be added and we obtain 

V> 2 = (1 - G tP)- 1 G tP(l - P 34 )^i = G fP(l - P 34 )Vi (3.23) 

with 

f = t l2 + t 12 PG f. (3.24) 

The above equation, Eq. (|3.24|) , defines the T-operator for the two two-particle subclusters, 
which interact only internally but not with each other. The desired second Yakubovsky 
equation is given by Eq. ( |3.23|) . For the sake of completeness we give also the final expression 
for the total scattering state 

*(+) = (1 + P - P 34 P + P)((l - P 34 )^ + </> 2 ) (3.25) 

In order to define an optical potential for the scattering of a nucleon from 3 He, 
we have to separate in Eq. ( |3.19|) the propagation in the 3 He subcluster channel from the 
propagation, where the three nucleons 1,2, and 3 are in intermediate scattering states. In 
order to achieve this, we have to reconsider the definition given in Eq. ( |3.17| ) in the following 

way: 

(1 - G ti 2 P)~ 1 G ti 2 P = (1 - GV 12 P)- l GV 12 P 

= (E-H -V 12 (1 + P))- 1 V 12 P 

= G TP. (3.26) 

Here we introduced the resolvent operator 

G= (E + ie-H -V 12 ). (3.27) 

In the four-body system the convenient 'odd man out' notation conventionally used in the 
three-body system is not applicable and what there appears 'natural' as choice of the ar- 
bitrarily singled out pair, namely 1 = (23), is not being applied to the four-body system. 
Instead, we have the pair (12) as the 'first' pair. Again, the specific choice of the 'first' pair 
is irrelevant and just a matter of convenience. 

By the very definition of the target state \<p) and its Faddeev component \<fi F ) given 
in Eq. (|3.15| ) one has 

(E He - /iq(123) - V 12 (l + P))\<P F >= (3.28) 

and 
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< <P\(E He - fc (123) - V 12 {1 + P)) = 0. 



(3.29) 



Here the 3N binding energy is given by Eu e , and ft-o(123) stands for the internal kinetic 
energy of the nucleons 1, 2 and 3. The above two equations indicate that there are different 
left and right eigenvectors, and thus we obtain 



(E + ie-Ho- "Ml + P))- 1 = |0 F > g — * < <j>\ + continuum , (3.30) 

< <p\Q > 



where 

9o = -TT-- k 2 ~^ (3-31) 

E + xe - E He - ^u 2 

is now the single particle propagator in the 3 He subcluster channel, while in Section II it 
was used for describing the propagation in the deuteron channel. Using Eq. ( |3.30| ) we now 
can cleanly separate this propagation from the rest 

G TP = \<P F > g - 1 < <p\V l2 P + G T C P. (3.32) 



Inserted into Eq. ( |3.19|) we arrive at the final expression for the Yakubovsky component fa: 

fa = $ F + |0 F > 9o— r^— < <f>\V 12 P\(-P 34 fa + fa) > +G T C P(-P u fa + fa). (3.33) 
< old) > 



Obviously, the second term on the right hand side provides an asymptotically purely outgoing 
wave carrying the elastic amplitude for p- 3 He scattering 

M = 1 < $| V 12 P(-P u fa +fa)>. (3.34) 
< <p\<i> > 



Here ($| = (0|(u| is the on-shell channel state. This form is clearly not the standard one and 
we sketch briefly the link to the more familiar expression. To achieve this, we first rewrite 
Eq. (pi as 

< 4>\4> F > M = < §\v l2 P(-Pz A fa + fa) > 

= <<5>\(E-H - V 12 )\(-P u fa + fa) > 

= < $\(E - H - V 12 )\(fa 2 - fa) > . (3.35) 

The second equality corresponds to Eq. (|3.29|) for ($| as introduced above and the 4N kinetic 
energy H . The third equality is due to Eq. ( |3.9| ). Inserting now the definitions of fa 2 and 
fa, Eqs. ( |3.3|) and ( |3.11|) , we obtain 

< <j)\(f) F > M = < -H - V l2 )G V 12 y {+) > - < ®\(E -H - V 12 )G tPG V 12 \¥ +) > 
= < ®\V 12 \¥ +) > - < <S>\V 12 G Q V 12 \¥ +) > - < <$>\V l2 PG Q V l2 \¥ +) > . (3.36) 



Now we use Eq. ( |3.1| ) and get together with the explicit definition of the permutation operator 
as given in Eq. (|2.14j ) 
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< > M =< <S>\V 12 G (V 14 + V 24 + V 34 )|^ (+) > . (3.37) 

Due to the antisymmetry of $ with respect to particles 1, 2 and 3 and of 
(V14 + V24 + V34) l 1 !'^^) this can be rewritten as 

< <p\<p F > m = 1/3 < $|(y la + v 23 + v 31 )G (Vi 4 + v 24 + v 34 )|^ (+) > 

= 1/3 < $|Vi 4 + V 24 + V 34 |^ {+) > . (3.38) 

In the last equality we took advantage of ($|(Vi 2 + V23 + V 3 i)G = ($|. Finally we note that 
= 1/3 for a normalized state This concludes our verification that the amplitude 
M given in Eq. (|3.34j) is indeed the desired elastic scattering amplitude for p- 3 He scattering. 

Starting from Eq. ( |3.33| ) and Eq. ( |3.23| ) one can already work out the optical po- 
tential formalism for p- 3 He scattering. We do not want to do this here but rather in addition 
separate off the propagation of two deuterons in Eq. ( |3.23| ). This allows to derive coupled 
equations for p- 3 He and d-d scattering. Very likely this will accelerate the convergence of the 
expansions of the resulting 'optical' transition potentials. In order to achieve the separation 
of the d-d channel we proceed analogously to the separation of the p- 3 He channel. First we 
note that in analogy to Eq. (3.26) 

G fP =(E-H - V 12 (l + P))- l V 12 P. (3.39) 

Now we need to introduce the analog to 'Faddeev' components for the d-d channel, which 
are somewhat unfamiliar. The two uncorrelated deuterons 0^(12) and 0d(34) obey the 
Schrodinger equation 

(/io(12) + Vi 2 - E d + /i (34) + V 34 - £ d )0 d (12)0 d (34) = 0, (3.40) 
where ho{ij) are internal kinetic energies only. Let us call 

4>dd = \4>d >12 \4>d >34 • (3-41) 
In view of Eq. ( ^.39| ) we can rewrite Eq. (|3.40| ) as 

< <p dd \(2E d - h (12) - h {M) - V 12 (l + P)) = 0. (3.42) 
The integral form of Eqs. ( |3.42p or ( |3.40| ) then reads 



>dd 



2E d -h {12)-h Q {34) 
The state <j) dd can now be decomposed in a Faddeev-like fashion as 



(Via + V 34 )<^ • (3.43) 



(Pdd = Vdd 
with 



/L2 + 0£f 4 . (3-44) 



& = & = 2E d -h (12)-h (34) Vl ^ dd (3 ' 45) 



and 
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= Hm- (3-46) 

It follows that 

(2E d - h (12) - ha(te) - V l2 (l + P))<p F dd = 0. (3.47) 
Using Eqs. ( |3.42| ) and (|3.47 ) we can separate the d-d channel from the rest 

G fP = \<f> M > g f 1 < ^ dd \ Vl2 p + Go f C p. ( 3.4 8) 

< <Pdd\<Pdd > 

Here the propagator in the d-d channel is explicitly given by 

1 

it- — 

2m 



where v is a suitable Jacobi momentum for the relative motion of the two deuterons. 
Insertion of Eq. (ggg ) into Eq. (|3~23|) gives for the second Yakubovsky equation 

^2 = \<f> M > gf —, ). F - < 0^12^(1 - P 34 )^i + G T C P(1 - P 34 )Vi- (3.50) 

The asymptotic form for the d-d channel follows in a similar fashion as for Eq. ( |3.33|) and 
the resulting transition amplitude into the d-d channel is given by 

M dd = - ] < $ dd \V 12 P(l - P 34 |^i >, (3.51) 

where & dd denotes the d-d channel state 

®dd = \<Pdd > |v > (3.52) 

Again, this is not the standard form for the transition amplitude. We want to briefly sketch 
the relation to the standard form. The basics steps are very similar to the ones leading to 
the elastic amplitude given in Eq. (|3.38|) . We start from 

< (j) dd \(j) M >M dd = < <& dd \V l2 P{l - P 34 |^i > 

= < <$> dd \{E -H - V 12 )(l - P 34 M > 
= < $ dd \(E -H - V 12 )\ip 12 -ip 2 > 

= < $ dd \(E -H - V 12 )G V 12 \¥ + ^ >-< $ dd \(E -H - V 12 )G tPG V 12 \¥^ > 

= < $dd\V 12 \V M > " < $dd\V 12 G V 12 \¥^ > - < $dd\V 12 PG V 12 \^ > 

= < ® dd \V 12 Go(V 13 + V u + V 23 + ^ 24 )|* (+) > 

= 1/2 < $ dd \(V 12 + V U )G (V 13 + V u + V 23 + V 24 )|^ (+) > 

= 1/2 < <S> dd \(V 13 + V u + V 23 + \/ 24 )|^ (+) > . (3.53) 

Since {<fidd\<j>dd) = 1/2 we confirm the standard form. 

Now we are ready to derive coupled equations for the p- 3 He and d-d channels 
starting from Eqs. (|3.33|) and (|3.50|) . As a first step we introduce coupled amplitudes for 



12 



general Jacobi momenta u and v corresponding to the driving terms either in the p- 3 He or 
the d-d channel, 

= ®u+ G T C P(-P 3 4^, U + r 2 ,u) 
r 2 , u = G T C P(1 - P 34 )^, u (3-54) 

and 

< V = G T C P(-P 34 < V + < V ) 

<v = + G f c P(l - P 34 )< v (3.55) 

Here $^ dv = 0^|v) is the 'Faddeev' component in the d-d channel together with the relative 
momentum eigenstate |v) of the two deuterons. The second set of equations, Eqs. ( p.55|) , 
is necessary to derive coupled channel equations. With the above definitions it follows from 
Eqs. and (^50|) 

+ / d 3 <flf(v) } < ^ d \V 12 P(l - P 34 )|Vi > (3-56) 

and 

V>2 = V2,uo + J d3u ^ J lu9oH < -^-p ~ < $u|^12P(-P 3 4^1 + ^2) > 

+ / ^ v flf(v) — ■ - * ~ < $f|Vi 2 P(l - P 3 4)|^l > • (3.57) 

•/ < <pdd\(pdd > 

Here we denote the initial on-shell relative momentum by Uo in contrast to the previous use. 
The right hand sides contain the transition amplitudes from the initial channel p- 3 He to the 
channels p- 3 He and d-d. Starting from these equations we can easily obtain the following 
set of coupled equations for half-shell transition amplitudes: 

< $u'|^12P|(-P34^1 +^ 2 ) > = < $u'|Vl 2 P|(-P34^, U0 + ^, U J > 
+ Jd 3 U < $ u ,|Vl 2 P|(-P 34 ^ iU + 4> C 2,u) > ^O(U) <0| 1 0F> < $u|V 12 P|(-P 3 4Vl + fo) > 

+ fd 3 v< $ u ^i 2 P|(-P 3 4< v + < v ) > ffoV) , ), F ~ < K d \Vi 2 P(l - PmM > (3.58) 
and 

< ^\V 12 P(1 - P 34 )|Vl > = < $V \Vl2P(l - P^Miuo > 

+ Jd 3 U< ¥J\V l2 P(l - P 34 M,u > go(u) < ^ > < $u|^l 2J P(-A4^1 + ^2) > 

+ /d 3 ^ < <|V 12 P(1 - P 34 )|< v > SoV) — -4rpr- < Oi 2 P(l - Pu)\1>i > (3.59) 
J < <Pdd\(Pdd > 

Defining transition potentials as 
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V u , u = < $u'|V 12 P|(-P3#l C ,u + V>2,u) > 

V v , v = < < \VuP(l ~ PuMtv > 

Vufy = < $u'|^l 2 P|(-P34< v + < v ) > 

K'u = < ® d J\V 12 P(l - PuWi,u) >, 
we can rewrite Eqs.( |3.58 ) and ( |3.5S ) in a more transparent form: 
<®u>\V 12 P\(-P 34 ip 1 + ij 2 )> = v u > y 



(3.60) 



+ 



'u'uo 

j3 , 



, ePuV u , u g {u) — p < ^n\V 12 P\(-PMipi + ip2) > 
J < (f>\(p p > 

+ J d 3 vV u < v i ^ ' 



7 ^<^|y la p(i-p 34 )|V> a > (3.6i; 



< |y 12 p(i 



^34)1^1 > = K',u 

+ / d 3 «V v / u # (u) 



< 



+ / f#vV9 (v 



< 



^ < $u|V 12 P|(-P 34 ^l+^2) > 



<< d |V 12 P(l-P 34 )|^i> (3.62) 



These are the coupled integral equations for the two-body fragmentation channels p- 3 He and 
d-d and constitute the main result in this study. The solutions lead directly to the elastic 
amplitude M and the transition amplitude into the d-d channel M^d as given in Eqs. (|3 . 34 ) 
and (|3T51| ) . 

The remaining task is to determine the amplitudes ip° and ipf as given in Eqs. (|3.54|) 
and (p.55| ). The idea and hope is that at high enough energies this set of amplitudes can be 
obtained by iteration and that only the first few terms will contribute. To be explicit let us 
consider the second order iteration of those equations, which is given by 



ft, 



2,u 



+ 



**\,( G T C P(-P 34 )<^ 
J + vG T c P(l - P 34 )$£ 
G T C P(-P M )G T C P(-P 34 )^ + G T C PG T C P(1 - P 34 )$£ 
G T C P(1 - P u )G T c P(-P m )$f 



(3.63) 



and 



ft, 



2,v 



+ 



(3.64) 



\ /G T c P$^ iV 

GoT c P(-P 34 )G r c P$^. v \ 
G T C P(1 - P 34 )G T C P$^ J 

Inserting this result into Eq. ( |3.60| ) we find for example the transition potential for p- 3 He to 
p- 3 He scattering in this order of approximation 

V u>iU = <® u ,\V 12 P\(-P 3A )\$F> 

+ < $ U ,\V 12 P(-P U )G Q T C P(-P U ) + V 12 PG f c P(l - P 34 )|$£ > 

+ < ®u>\V 12 P(-P 3i )G T c P(-P u )G T c P(-P u ) + y 12 P(-P 34 )GoT c PGoT c P(l - P 34 ) 

+ y 12 PG T c P(l - P 34 )G T C P(-P 34 )|$^ > . (3.65) 
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The other transition potentials can be obtained in a similar fashion and are given in Ap- 
pendix B. 

In order to evaluate the transition potentials we have to determine T c and T c . 
According to the definitions given in Eqs. (|3.32| ) and (|3.48|) as well as the defining equations 
( p,18| ) and ( |3.24| ) for T and T, we obtain after some algebra 



G T C 



> 



1 



< 



> 



< 



Gotu + GoT c PGot\2 



(3.66) 



and 



GnT c 



1 - 



< (f>dd\ 



j dd 



> 



< 4>dd\ 



Goti2 + GoT c PGoti2. 



(3.67) 



As a consequence of the uniqueness of these equations and also from their very definition 
(the latter one requiring some algebraic manipulations) one can conclude that 



< 4>\G,T C 







and 



< (f> dd \G T c = 



(3.68) 



(3.69) 



This excludes intermediate propagation of 3 He or d-d states in the amplitudes given in 
Eqs. (|3.63|) and ( p.64j ). Of course, this should be the case. From solving 3N Faddeev 
equations at intermediate energies [17], we know that the multiple scattering series converges, 
except for the 3 He quantum numbers J 71 " = l/2 + . This divergence, however, is due to the very 
existence of 3 He. Due to the special driving terms in Eqs. ( |3.66| ) and ( |3.67| ) this divergence 
is removed. Therefore, there is a good reason to assume that Eqs. ( |3.66|) and (|3.67| ) can be 
successfully iterated and that very low orders are sufficient. With the abbreviations 



A = 1 



> 



< 



> 



< 



and 



A. 



dd 



1 



J dd 



> 



1 



< (f>dd\ 



'dd 



> 



< 4>dd\ 



we consequently approximate GqT c and GqT c to second order in tu as 

GqT c — > AGoti2 + AGo^i2-PGoti2 

and 

GqT c —* AddGotu + AddGotuPGotu. 



(3.70) 



(3.71) 



(3.72) 



(3.73) 



When inserting these expressions into Eq. ( |3.65| ) we obtain for the transition potential V u i 
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V u , u = <<f> u ,\V 12 P{-P 34 )\^> 

+ < $ u ,\V 12 PP 34 AG t 12 PP 34 + V l2 PA dd G t 12 P[l - P 34 )|$u > 
+ < ®u>\V 12 PP u AG t 12 PGot 12 PPu + V 12 PA dd G t 12 PG t 12 P(l - P 34 )|$u > 
- < $ u ,\V 12 PP 34 AG t 12 PP 34 AG t 12 PP 34 + V 12 PP 34 AG t 12 PA dd G t 12 P(l - P 34 ) 
+ V 12 PA dd G t 12 P(l - P 34 )AG t 12 PP 34 \^ > . (3.74) 

In the parts of V u ' u which are second order in T c and T c we kept only the first order parts 
of Eqs. ( |3.72| ) and ( |3.73j ) which are linear in ti 2 . 



As one of several possible simplifications in the above expression, we note that P(l — P 34 ) = 
(1 — P 42 )P and (1 — Pi 2 ) yields a factor of 2 upon its application to t 12 . Thus we obtain for 
the transition potential 

K'u = <$u'|V 12 P(-P 3 4)|^> 

+ < $ u ,|^ 12 PP 34 AG t 12 PP 3 4 + 2V 12 PA M G t 12 P\^ > 

+ < ®u>\V 12 PP 34 AG t 12 PG t 12 PP 34 + 2V 12 PA dd G t 12 PG t 12 PK > 

- < $ u <\V 12 PP 34 AG t 12 PP 34 AG t 12 PP 34 + 2V 12 PP 34 AG t 12 PA dd G t 12 P 

+2V l2 PA dd G Q t l2 PAG Q t l2 PP 34 \<b F u > . (3.75) 

The evaluation of this expression requires the know-how of handling all the separate sub- 
cluster problems: The 3 He bound state and its Faddeev components, the deuteron and the 
Faddeev component to the two deuteron subcluster problem, the NN t-matrices and above 
all the (3+1) and (2+2) subcluster kernels G t 42 P and G ti 2 P. In addition there is the 
P 34 permutation operator to be considered. All those mathematical structures have been 
successfully handled in previous work |i6||T7|l , and thus the evaluation of this transition 
potential should pose no additional difficulty, although it is certainly a computational chal- 
lenge. Similar expressions can be derived for the remaining transition potentials. They are 
given in Appendix B. 

A close inspection of Eq. Q3.75| ) shows that most of the terms will likely not 
contribute at higher nucleon projectile energies, since they are exchange terms, where the 
projectile momentum probes bound state wave functions in a region where they are already 
very small. But there are of course also 'tp'- type structures, like part of the fourth term 
in Eq. ( ^.75| ), which will dominate. We leave a more detailed study of Eq. ( |3.75| ) to future 
work. The importance of different terms should become apparent when being supported by 
numerical realization. Nevertheless we would like to emphasize that the presented approach 
is systematic. It relies on a mathematically and physically well founded basis. Scattering and 
bound state structures are treated on the same footing, and antisymmetrization is included 
fully. The internal validity of our approach can then be checked numerically by adding 
further terms of the expansions. There will be no free parameters in the calculation, once a 
certain NN force has been selected. 

At a later stage it might be also of interest to study the properties of those tran- 
sitions potentials with respect to their spin-dependencies and locality versus nonlocality. It 
would be a surprise if they would have much in common with the often used phenomeno- 
logical Wood-Saxon type expressions. 



16 



IV. THREE- AND FOUR-BODY FRAGMENTATIONS IN PROTON- 3 HE 

SCATTERING 



As alternative to an algebraic derivation as pursued in the previous section, we 
would like to arrive at the set of coupled integral equations for the breakup process as given 
in the Yakubovsky scheme starting from a graphical approach. The resulting equations are 
exact and can be derived rigorously in an algebraic manner. 

The complete breakup process initiated for example by nucleon number 4 striking 
a 3 He target composed of nucleons 1,2, and 3 is given by the infinite sequence of processes 
depicted in Fig. I. Here to the very right of each diagram represents the 3 He target 
ground state and the dashed lines stand for NN interactions. Each diagram has to be read 
from right to left and each has to start with an interaction between the projectile and one 
of the constituents of the target. After this initial interaction of the projectile, arbitrary 
interactions between all four particles have to occur. Clearly, this comprises all possible 
interactions and intermediate free propagations of the four nucleons. The superscript on the 
breakup operator indicates that here particle 4 has been singled out as projectile. In 
order to achieve full antisymmetrization, only the interchange between the projectile and the 
single target nucleons have to be considered, since the target ground state is assumed to be 
already antisymmetrized. Thus, the fully antisymmetrized breakup operator U is given by 
the set of diagrams in Fig. 2. There, the terms representing nucleons 1, 2 or 3 as projectiles 
enter with a negative sign, since only one transposition is necessary to interchange nucleon 
4 with one of them. 

(k) 

For each of the three terms Uq , k — 1, 2, 3 in Fig. 2 expansions corresponding to 
Uq^ as given in Fig. 1 can be written down. The fully antisymmetrized breakup operator 
Uq can be decomposed into 6 Faddeev components according to the last pair interaction on 
the left, 

Uo = Y.U?y (4.1) 

i<j 

By inspection of the various Born series represented in Fig. 2, one can read off the different 
Faddeev components. For example, the component U® 4 is given by amplitude 

C/ 3 °4 = ^34(<f 4 - $ 3 ) + V U Gq E U?j- ( 4 - 2 ) 

i<j 

Here |<J> 4 ) = |0(123))|u) 4 and |$ 3 ) = P 34 |<E> 4 ). As usual one can sum up V 34 to infinite order 
into t 34 : 

*<i,*i^34 

= t 34 ($ 4 - $ 3 ) + t u G (U° 4 + C/ 2 ° 3 ) 

+i 3 4Co(C/ 1 4 + E/? 3 ) + t u G U° 12 . (4-3) 

In the second equality we already group the terms in a suggestive manner, namely such that 
the pair indices in the last three terms define two different three-body subclusters, namely 
234 and 134, and one 2+2 fragmentation 12-34. 
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From here on we use the identity of the particles, which simplifies matters consid- 
erably and allows to restrict the discussion always to one type of amplitude, either from the 
partition 3+1 or from the partition 2+2. All remaining amplitudes are obtained by suitable 
permutations of the particles. Like in the previous section we define a (3+1) Yakubovsky 
component 

[7° = t u G (U° 4 + f/ 2 ° 3 ) (4.4) 

and a (2+2) Yakubovsky component 

U\ ee t u G U° 12 (4.5) 

as two independent amplitudes. Using the identity of the particles, the (3+1) component 
becomes 



Ui — t34Go(— -P23 — F > 2i)U [ 



34 



= t u G (-P 23 - P 24 ) [t 34 ($4 - $3) + Ul - P12U? + U°] . (4.6) 
It is convenient to introduce a permutation operator for particles 2,3, and 4 

-^234 = -^23-^34 + -^24-^34; (4-7) 



since 



Then 



P uU ° M = -U» A . (4.8) 



Ui — ^34^0^234^34(1 — P 3 a)$A 

+tuG G P 2 u{l - P ia )Ui + t 3 ,G P2 3 4Ul (4.9) 



Analogous to Eq. (|3.26|) we sum up all interactions in the 234 subsystem to infinite order 
and obtain 

U°i = TG P234*34(1 - P 34 )$4 + TGoP2u(-PnUl + C/ 2 °). (4.10) 
The above equation defines the (off-shell) three-particle operator T as 

(1 — £34^0^234) 1 ^34G r o-P234 = TG0P234, (4-11) 

which is equivalent to the integral equation 

T = t u + t u G P 2 uT. (4.12) 

For this off-shell three-body operator T the same notation is used as in Eq. ( |3.18| ), since it 
represents exactly the same quantity, however it is expressed in different particle numbers. 
Similar steps as those described above lead to 

U° = TG Pt u (l - P 34 )$4 + TG P(1 - P U )U°, (4.13) 



18 



where P is given in Eq. ( |3.20|) and 

T = £34 + t^G^PT. 



(4.14) 



Again the same remark as above applies concerning T. The coupled set of Eqs. ( [4.10| ) and 
Q4.131 ) for Ui and U® are the Yakubovsky equations containing (3+1) and (2+2) subcluster 
T-operators. Once solved, the Yakubovsky components and U® are sufficient to generate 
the full breakup operator U . 

Summing all Faddeev components as given in Eq. (|4.1j ) and taking advantage of 
the identity of the particles yields for the breakup operator 



u = (i + p 234 + P- p^aWI^ 



with U® 4 given in Eq. ([Oj] 



Uu = *34(1 - A 4 )$4 + (1 - PM + Ul 



(4.15) 



(4.16) 



At intermediate energies one can expect that the lowest order terms in the NN t-matrix 
should be sufficient. When inserting the driving terms of Eqs. ( |4.12 ) and Q4.14 ) into the 
Yakubovsky equations ( |4.10 ) and ( 4.13Q , we obtain in lowest order 



U\ ~ ^34^0^234^34(1 — Pm)^a 
C/ 2 ° « t M G Q Ph A {l - P 34 )$4. 

Thus, in lowest order the breakup operator Uq is given as 



(4.17) 



U = (1 + P 23 4 + P - 
*34(1 - P 3 4)$4 



-P12-P234) 

f (1 " Pl 2 )t34G P234t34(l - P 3 4)$4 + ^34^0^34(1 - ^34)^4 



(4.18) 



This expression is manifestly antisymmetric in all four particles, since the square bracket is 
separately antisymmetric in the pairs 12 and 34, as is obvious for the first and second term. 
Similarly, this can be seen for the third term by noting that ^£34 (1 — -P34) = £12(1 — Pn)P- 
Now P applied onto <J> 4 yields — $ 2 , which is antisymmetric in 34. The total antisymmetry 
can then be checked using the permutation operators in front of the square bracket. 

As an example for the actual evaluation of the terms in Eq. ( |4.18j ) let us consider 
the term of first order in the NN t-matrix. It is a simple exercise to show that 



1 + P 234 + P- Pl 2 P 2 34)t34(l - P 34 )$4 

= £ 14 (1 - Pl 4 )$4 + £24(1 " P 24 )$4 + £34(1 - P 34 )$4 

-£ 12 ((1 - P 12 )$ 2 - £ 13 (1 - Pl 3 )$l - £ 23 ((1 - P 23 )$2- 



(4.19) 



The breakup amplitude results by acting from the left with momentum eigenstates for four 
free particles. They can be represented by 3 Jacobi momenta suitably chosen according to 
the pair interacting in the various NN t-matrices. In addition, the various channel states 



by choosing 



1,4 require adequate choices of Jacobi momenta. For instance $4 is best described 
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P = ^( k l - k 2) 

q=^(k 3 -^(ki + k 2 ) 

r = ^(k4-i(k 1 + k 2 + k 3 )). (4.20) 

Then 

< pqr|$ 4 >= 5(r-r o )0(pq), (4.21) 

where ro is the initial projectile momentum (in the zero total momentum frame) and 0(pq) 
stands for the 3N bound state. For the sake of simplicity we drop spin and isospin quantum 
numbers. In order to describe for instance £34 it is convenient to introduce a second set of 
Jacobi momenta, which singles out the relative momentum among particles 3 and 4: 



Pi = 


i(k3 


-k 4 ) 


qi = 




-^(k 3 + k 4 ) 


ri = 


4 (kl 


-i(k 2 + k 3 + k 4 ) 



(4.22) 



Then one easily derives 



2 1 ^ 

< Piqiri|pqr > =b(p- -r a + -qi) x 



3 2 

2 2 4 1 1 

<Kq - 3PI + gQl + g r l)^( r + Pi + 2*11 + 3 r l) 

2 1 2 2 4 

*(pi + g 1 " - 2-q)<Kqi + 3P + g-q + g r ) x 

5( ri -p + iq+ir). (4.23) 



Furthermore, one has 



3 2 

< PiqiriMpiqirl >= t^Pi, pi; E - —q\ - — r?)5(qi - ql)<*(ri - r[). (4.24) 
Therefore, using Eqs. ( |4.21| ) and ( |4.24[ ) one easily finds 

< piqir 1 |t 34 (l - P 3 4)$4 > = 

1 1 3 2 2 2x /2 1 2 2 \ . 

2 t 34 (pi, -r - -qi - -n ; E - — qi - — rj* - - qi , --r - * - -n j . (4.25) 

This term is maximal under the condition, that both arguments in are zero. This leads to 

4 

qi = -gro 

n = -~r (4.26) 
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Noting that tq = f k^, where k; a fe is the projectile momentum in the laboratory system, 



one easily deduces by using Eq. Q4.22|) , that the conditions Eq. ( |4.26|) are equivalent to 
ki = k2 = 0. The condition of zero momenta for the spectator nucleons is usually called 
quasifree scattering. 

Next we consider the arguments of the NN t-matrix £ 34 in Eq. ( |4.25| ). Using 

E = E He + Ar 2 = d + ± q f + ±rl (4.27) 
6m m Am 6m 



where Ejj e is the binding energy of 3 He, one finds under the condition given in Eq. Q4.26 ) 



*34(pi, -^r ; E He + -(^r ) 2 ). (4.28) 
6 m 6 

Therefore £34 is on the energy shell except for the negative binding correction En e of the 
3 He target. 

Going away from the quasielastic peak one can probe the off-shell NN t-matrix 
and the 3 He target wave function. At the same time however, the other terms in Eq. (|4.19|) 



will also contribute. Each term in Eq. ( |4.19|) expressed in the adequate Jacobi momenta will 
have a form similar to Eq. (|4.25|) . Corresponding spectator momenta set to zero have the 
effect, that one term will dominate. The other ones are suppressed, since the arguments of 
4> then differ from zero and the t-matrices are off-shell. Numerical studies are required to 
learn about interferences among the six terms if one is away from the peak conditions. One 
can expect to extract interesting information on the ground state wave function <fi of 3 He 
and the off-shell NN t-matrix. Both are nowadays readily accessible for realistic NN forces 
and can therefore be tested. 

The next step is to evaluate the terms second order in t 34 as given in Eq. (|4.18|). 



This is more complicated, since it involves the free propagator, which together with the 
permutation operators leads to logarithmic singularities. However, their treatment is known 
from the three-body system |l7j, and thus these terms can also be numerically determined. 



We expect that with increasing energy the multiple scattering series should terminate quickly 
and therefore the reaction as well as properties of 3 He can be tested systematically and in a 
rigorous way. This will be left to future work. 

A final step is the formulation of the transition to the three-body fragmentation 
channels d+p+p. Again we start by choosing particle 4 as projectile and display the first 
terms of the infinite series of processes in Fig. 3. Here one has to note that the last interaction 
to the left cannot take place between particles 1 and 2, since V12 is already taken into account 
in the deuteron ground state wave function <p d {12). We recognize the same series of processes 
as in Fig. 1, except that the last interaction V12 is excluded. The antisymmetrization in the 
initial state leads to the symmetrized transition operator 



U d =uf-Y,U$\ (4.29) 

i=l 



where refers to the projectile i. For 2 = 1,2 there occurs just one term with V12, namely 
Vi2</>i, i = 1,2. Thus we have 



u i = -V i3 (<l> 1 + fa) + J] u ir ( 4 - 3 °) 

i<j,ij^l2 
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For the on-shell process the potential V12 can be replaced by Go 1 - Finally, permutation 
operators can be introduced leading to 

U d = -G \P 1A + P 24 )$ 4 + (1 + P - Pi2P)Ul 4 . (4.31) 

Here is given in terms of the two Yakubovsky amplitudes and U® of Eq. (|4.16| ). Again 
the lowest order terms in Eq. ( |4.17| ) are expected to be sufficient at intermediate energies. 
This remains to be verified in a numerical study. 



V. SUMMARY AND OUTLOOK 

Elastic proton-deuteron scattering is formulated within the framework of Faddeev 
equations to achieve the form of a one-body equation with an optical potential as driving 
term. The exact equation defining this potential is assumed to be solvable at intermediate 
energies by a low order expansion in the NN t-matrix. This remains to be verified numer- 
ically. The presented formulation is a systematic expansion in two-nucleon t-matrices and 
treats the antisymmetry among all three particles exactly. A standard approach to derive 
microscopic optical potentials at intermediate energies is the 'spectator expansion' of mul- 
tiple scattering theory. We apply this method in the case of p-deuteron scattering in order 
to illustrate similarities and differences to the treatment within the Faddeev scheme. 

In the four-nucleon system we formulate an exact set of coupled equations for p- 
3 He and d-d scattering within the framework of the Yakubovsky equations. The transition 
potentials between those two channels are again approximated at intermediate energies in a 
low order expansion in the NN t-matrix. This expansion is systematic and antisymmetriza- 
tion among all four particles is treated exactly. The numerical realization as well as the 
internal check of convergence with respect to higher order terms should provide interesting 
insight into the reaction mechanism and should lay a firm ground to test various physical 
assumptions, like the application of free NN forces in such a reaction as well as the properties 
of the 3 He wave function as resulting from solving the Faddeev equations based again on 
NN forces. 

We also derived exact equations for the p- 3 He induced three- and four-body frag- 
mentation processes. Here we again concentrate on the lowest order terms in a NN t-matrix 
expansion, which we expect to be valid at intermediate energies. The present work lies a for- 
mal ground for numerical investigations, which are planned. The expectation is that at least 
at intermediate energies the four-nucleon scattering problem can be numerically controlled 
in a reliable manner, using the outlaid systematic t-matrix expansion. At lower energies such 



an expansion fails, as is already known from the 3N system |17| . In the low energy regime 
an exact solution of the four-nucleon Yakubovsky equations appears not to be feasible at the 
present time. Since the proposed scheme does not contain any adjustable parameters once a 
realistic NN force has been chosen, its numerical realization and comparison to experiment 
should provide interesting insight into our understanding of such reaction processes and the 
properties of the 3 He wave function. 
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APPENDIX A: THE SPECTATOR EXPANSION FOR P-D SCATTERING 

In this Appendix, the most simple system for nucleon-nucleus scattering, the p-d 
system, is used to explain the steps involved in the spectator expansion of multiple scattering 
theory, on which one model for the nucleon-nucleus optical potential is based ||. Since the 
p-d system is the simplest system for nucleon-nucleus scattering all steps can be clearly 
carried out and the result can be compared to the exact Faddeev framework of Section 
II. We especially want to emphasize the treatment of the antisymmetrization and show 
the inherent limitation in spectator expansion. For distinguishable particles the spectator 
expansion is carried out for three particles in Ref. |]20| . 



In the general derivation of the optical potential for scattering of a nucleon from 
a composite nucleus, the protons and neutrons are treated as distinguishable particles. In 
this Appendix we will follow this general practice. Numbering the two protons as particles 
1 and 2 and the neutron as particle 3, the scattering state initiated by proton 1 scattering 
from a deuteron composed of nucleons 2 and 3 is given by 

where G the full resolvent operator and 

|$i >= |fo(23) > |q„ >i • (A2) 

The scattering state antisymmetrized in the two protons is then explicitly given by 

*(+) = (1 - P 12 )^ +) (A3) 

and enters the matrix element for elastic p-d scattering as 

M =< $ q ,|(V 12 + y 13 )|^ (+) >=< $ q /|(Vi2 + V 13 )(l - P 12 )\v[ +) > . (A4) 

Here |3> q ') is defined analogously to Attempts to derive Lippmann-Schwinger equations 
for the fully antisymmetrized scattering state given in Eq. ( |A3[) did not lead to expressions, 
which are applicable in practice EHf . 



As an aside, we would like to mention a fact, which may not generally be known. The 
fully antisymmetrized state \E ,< - + - ) obeys the same Lippmann Schwinger equation as 
namely 

*W = + Gi(Vi3 + V 12 )y{ +) . (A5) 

The simple reason is that due to Lippmann identities [p3|] , = Pi2^t obeys the homo- 
geneous version of Eq. (|A5|) . The same is true for a target composed of a general number of 
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A nucleons (N neutron and Z protons). This simply reflects the fact that Eq. ( |A5|) does not 
uniquely define the scattering state. Using Eq. (|A5|) for the state is only meaningful, 
if in the course of the solution antisymmetrization is always imposed, which is not easily 
implemented. 

Here we do not want to follow these thoughts. Instead we apply the operator P\ 2 in 
Eq. (lAl) to the left. This results in 



m = < $ q ,|(v 12 + v 13 )\v[ +) >-< q ,|(y 12 + v 23 )\^[ +) > 

= < $ q ,|(l - P 12 )V 12 \^ > + < $ q ,|Vi3lM +) > - < %>\Pi2V 23 \v{ +) > . (A6) 

The generalization to a number of A target particles (N neutrons and Z protons) can be 
easily obtained in a similar fashion. 

The first term in the second equality of Eq. ( |A~6|) imposes antisymmetry among 
the projectile target protons by antisymmetrizing the final channel state. The second term 
includes the proton-neutron interaction, which requires no antisymmetrization. The last 
term is an exchange term different from the first one, since the interaction V 23 acts inside 
the initial target state. As a note, this last term is always considered of higher order in 
applications of the spectator expansion model ||, and is thus neglected when only the first 
order term is considered. 

Since the physical interaction V 12 between the two protons has to be represented 
under all circumstances in physically allowed states, i.e. states which are antisymmetric 
under exchanges of (12), the operation (1 — P\ 2 ) provides just a factor of 2. Neglecting the 
last term in Eq. ( |A6| ) we are left with 

M -> 2 < $ q '|yi 2 |^i +) > + < $ q '|^i 3 |*i +) > • (A7) 

In the following we sketch the steps usually carried out in the spectator expansion. 
The Green's operator G\ in Eq. (|A5|) is projected onto a part propagating in the deuteron 
target state only, G\, and a remainder 

^(+) = + Gl(V 13 + V 12 )v[ +) + G{(V 13 + V 12 )^\ (A8) 

By introducing 

Xq +) = ^ + Gl(V 13 + V 12 )x^\ (A9) 
one finds together with G\ — J dV|$ q / > po(q') < $q| 

= X { £> + / d'q'x^goti) < d> q ,|(V 13 + V 12 )|*! +) > . (AlO) 



In view of Eq. ([A 71) it appears to be inevitable to treat the actions of V\ 2 and V13 separately 
and work with two coupled equations referring to the proton and neutron parts of M, re- 
spectively. However, doing so does not result in a single optical potential, which can serve 
as driving term for a one-body Lippmann-Schwinger equation for the elastic amplitude. 

From Eq. (|A10|) one can derive two coupled equations for the matrix elements 
involving V\ 2 and V13: 
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+ J d 3 q' < %\V 12 \ X ( ^ > po(q') < QjKVu + V 13 )\v{ +) > 
<$ q |V 13 |^ +) > =<%\V 13 \x^> 

+ J d 3 q' < ^\V 13 \x^ ] > </o(q') < %>\(V 12 + V 13 )\^ ] > (All) 

The next step is to evaluate Eq. (|A9|) . If one defines 

y liXq +) = r<|$ q >, (A12) 

where z = 1,2, one obtains for Eq. ( |A9| ) 

T, t % = Vu®* + VuG\ £ (A13) 

i=2,3 

Collecting all Tj$ q 's on the left side and inverting the result yields 

T < $ q = r i $ q + r i G5T i $ q (j ^ *), (A14) 
where the operators 7$ obey 

r< = V H + ViiGlTi . (A15) 

This result corresponds to the familiar Watson expansion Q. In first order one neglects the 
second part in Eq. (|A14|) , which describes the consecutive scattering of the projectile proton 
on both target nucleons. In this first order approximation, the only one which has been 
applied in practice, one has 

Ti% -> Ti %, (A16) 
and the optical potentials for pp and np scattering are given by 

< $ q |K 2 |x q t } > -> < $ q |r 2 |$ q , > (A17) 

< ^qlKalXqt } > -»• < $ qN$q' > • 

Denoting the two transition amplitudes in Eq. (|A11|) for simplicity by Z 2 = (^ql^lXq"]^) an d 
Z 3 = (^qjVislXq^) an d using a matrix notation, indicated by the 'bar '-operators, Eq. (|A11|) 
takes the form 

Z 2 = f 2 + f 2 g (Z 2 + Z 3 ) 

Z 3 = f 3 + ? 3 g (Z 2 + Z 3 ). (A18) 

It remains to solve Eq. (|A15 ), which takes the form 

n = V u + VuGrn - VuGln . (A19) 

The propagator G± still contains the interaction of the struck target nucleon i with the other 
target nucleon. The simplest way to solve Eq. ( |A19| ) would be to replace 



25 



Gi -> G • (A20) 

As an aside, expanding Gi = Go + G0V23G1 would lead to 'medium corrections' given by 
the potential V23. However, for the clarity of presentation, we do not want to pursue this 
further. In the approximation introduced by Eq. (|A20|) and introducing the free NN t-matrix 
tu Eq. ( |A19| ) simplifies to 

n = tu - hiP\n . (A21) 
Sandwiched between channel states this reads in matrix notation 

fi = hi ~ hmhi = (1 + hi9o)~%i ■ (A22) 

If one keeps the interaction V23 in G\ one faces a three-body problem, which 
has never been solved correctly in that context. The work in [^] provides an approximate 
treatment of the interaction of the struck target nucleon with the remainder nucleons treated 
through a mean field interaction. With the three-body technology available today, a correct 
treatment appears to be feasible and worthwhile. 

Continuing with the approximation of Eq. (|A20| ) and inserting the expression of 
Eq. (|A22|) into Eqs. flXIgp one finds 

Z 2 = in + h29oZ 3 (A23) 

Z3 = h3 + ^13^0^2- 

After a partial wave decomposition, these are simple one-dimensional integral equations. 
Their on-shell solution provides according to Eq. ( |A7| ) the physical transition amplitude 

M = 2Z 2 + Z 3 . (A24) 

The optical potentials occurring in Eqs. ( |A18| ) 



1-2,3 =< %\r 2l3 \%' > (A25) 

are the full-folding expressions of the NN t-matrices modified according to Eq. ( |A2 1| ) and 
integrated over the single nucleon density matrix generated from the deuteron wave function. 
The NN t-matrix contains the free 3N propagator Go = (E + is — 2- — -^q 2 ) and thus in 

2 

addition to the kinetic energy of relative motion within the pair, E -, also the kinetic energy 
of the pair as a whole in relation to the third particle, -^q 2 - In Ref. [0 the kinetic energy 
of the pair is treated within the full-folding model, while in Refs. |],[| this dependence is 
frozen into a constant. At this point we do not want to discuss the rationale for the two 
different treatments, but rather refer to the literature. 

It is appealing to compare the expressions of Eq. ( |2.27| ) for the optical potentials 
in the Faddeev treatment to the one of Eq. ( |A22j ) in the spectator expansion. However, this 



is quite hampered. In Eq. ( 2.27 ) a NN t-matrix t c occurs which is modified by a deuteron 
state for the same pair of nucleons as for the t-matrix, whereas in Eq. ( |A22|) the deuteron is 
composed by a different pair of nucleons. Furthermore, in Eq. ( |2.27|) there is the correction 
term (1 — Go^cT 1 ) and the factor 2 resulting from (1 — P\j). In the spectator expansion 
scheme the factor of 2 enters only in the final form of Eq. (|A24 ) and not in the integral 



equations (|A18 ), which are set up for distinguishable particles. 



In view of all these differences a numerical study would be of great interest to 
clarify at least in the three-body context the validity of the spectator expansion against the 
systematic Faddeev approach. 
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APPENDIX B: THE TRANSITION POTENTIALS 



The transition potentials entering the coupled integral equations for the p- 3 He 
and d-d channels are defined in Eq. ( |3.60| ). After considering the amplitudes ipf and ipf up 
to second order in T c and T c we find for the optical potential for elastic p- 3 He scattering, 
V u i u , the expression given in Eq. ( |3.65| ). In a similar fashion we obtain for the remaining 
transition potentials the following: 
The d-d to d-d transition potential is given by 

V v , v = < $$\V 12 P{1 - P U )G T C P\^ V > 

+ < ^\V 12 P(1 - P U )G T C P(-P U )G T C P\^ V > (Bl) 

The d-d to p- 3 He transition potential is given by 

Vtfv = < $ u '|Vi 2 PK )V > + < ^\V 12 P(-P U )G T C P\<I>^ V > 
+ < $u>\V 12 P(-P u )G T c P(-P u )GoT c P\^ v > 

+ < $ u >\V 12 PG f c P(l - P U )G T C P\^ V > (B2) 

The p- 3 He to d-d transition potential is given by 

K'u = < *$\Vi 2 P{l - P 34 )K > + < ^1^(1 " P3±)G T c P(-P u )\$f > 
+ < $^|V 12 P(1 - P M )G T C P(-P U )G T C P(-P U )\^ > 

+ < $^|V 12 P(1 - P u )G T c PG f c P(l - P 34 )K > (B3) 

After introducing for GqT c and GqT c the approximations given in Eqs. ( |3.72j ) and ( |3.73| ) 
as well as considering only terms up to second order in t\ 2 , we obtain for the transition 
potentials: 

K'v = < ^\V 12 P(1 - P u )AG t 12 P + V 12 P(1 - P 34 )AG t 12 PG t 12 P\^ d , v > 

+ < ^\V 12 P(1 - P u )AG t 12 P(-P u )AG t 12 P\^ dv > . (B4) 



and 



Kx'v = < $ u >\v 12 p\^ v > 

- < $u>\V 12 PP M AGot 12 P + V 12 PP M AG t 12 PG ti 2 P\^ dtV > 
+ < ^ u '\V 12 PP 34 AG t 12 PP u AG t 12 P\^ v > 

+ < ^\V 12 PA M G Q t 12 P(l - P u )AG t 12 P\^ v > (B5) 

V v , u = <$ d v 1\V 12 P(l-P u )\^> 

- < $^|y 12 P(l - P 3 4)AG ti 2 [1 + PG t 12 ] PP34K > 

+ < ^\V l2 P(l - P 34 )AGotl 2 PP34AG'otl 2 PP34| < ^ > 

+ < ^\V l2 P(l - P u )AG t 12 PA dd G t 12 P(l - P 34 )|$£ > . (B6) 



After the same simplifications leading to Eq. (|3.75|) we obtain the following expressions for 
the transition potentials: 
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VW = < ®P\2V 12 PAG t 12 P\<S>^ v > 

+ < ¥J\2V 12 PAG t 12 PG t 12 P\^ v > 

- < ¥ v d \2V 12 PAG t 12 PP u AG t 12 P\^ v >, 

V u , v = < $ u ,|y 12 P|$J d)V > - < $ u ,\V 12 PP u AG t 12 P\^ 

- < $ u >\V 12 PP 34 AG t 12 P[l - P u A]G t 12 P\^ v > 
+ < ^ u/ \2V 12 PA dd G Q t 12 PAG t 12 P\^ v > 

K'u = < ^\2V 12 P\$f >-< <S> d J\2V 12 PAG t 12 PP u \$F 

- < $ d J\2V 12 PAG t 12 P[l - P u A}G t 12 PP u \$F > 
+ < $ d v d \4V 12 PAG t 12 PA dd G t 12 P\^ > . 
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FIGURES 



FIG. 1. Born series for the four-nucleon breakup process p- 3 He — > pppn. Particle 4 is singled 
out as projectile. For further explanation see text. 

FIG. 2. The fully antisymmetrized four-nucleon breakup operator Uq. Since the target state is 
assumed to be antisymmetric only the exchange terms with the projectile have to be considered. 
This leads to four terms with each nucleon being the projectile. 

FIG. 3. Born series for the three-body fragmentation process p- 3 He — > ppd. Particle 4 is singled 
out as projectile. For further explanation see text. 
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